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Abstract. The energy input into the lower solar corona by flare evaporation events has been modeled according 
to the available observations for quiet regions. The question is addressed whether such heating events can provide 
the observed average level of the coronal emission measure and thus of the observed flux of extreme ultraviolet 
(EUV) and X-ray emission without contradicting the observed average power spectrum of the emission measure, 
the typical emission measure variations observed for individual pixels and the observed flare energy distribution. 
As the assumed flare height influences the derived flare energy, the mathematical foundations of nanoflare dis- 
tributions and their conversion to different height assumptions are studied flrst. This also allows a comparison 
with various published energy distributions differing in height assumptions and to relate the observations to the 
input parameters of the heating model. An analytic evaluation of the power spectrum yields the relationship be- 
tween the average time profile of nanoflares (or microflares), assumed to be self-similar in energy, and the power 
spectrum. We find that the power spectrum is very sensitive to the chosen time profile of the flares. Models are 
found by numerical simulation that fit all available observations. They are not unique but severely constrained. 
We concentrate on a model with a flare height proportional to the square root of the flare area. The existence of 
a fitting model demonstrates that nanoflare heating of the corona is a viable and attractive mechanism. 
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1. Introduction 

Many mechanisms have been proposed in the 
past for heating the solar corona (e.g. reviews in 
lUlmschneider. Rosner. &: Priesd. Il99l|) . Most popular at 
present are models based on either the release of magnetic 
energy or the dissipation of waves. Both scenarios were 
proposed soon after the discovery of the high coronal 
temperature by Edlen and Grotrian in the late 1930s. 
As the magnetic energy dominates the coronal energy 
density, the release of free magnetic energy, possibly 
present in the corona in the form of electric currents, 
is a very suggestive idea. However, it was soon realized 
that the electric resistivity caused by collisions between 
electrons and ions is extremely low in the thin corona, 
making this mechanism insufficient. Therefore, several 
wave heating mechanisms have been proposed. Most 
recent theories involve Alfven waves which can explain 
both the preferred heating of regions with high magnetic 
fields and solar wind acceleration (e.g. iMarsh fc Tul . 
IT997I). 

iGoldl l)l964f) seems to have been the first to suggest that 
the building up of free magnetic energy in the corona must 
be dissipated by flares when the resistivity becomes finite 
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by some inst a bility. The ide a was further develo p ed by 
iLevind (Il974tl. IParkerl lll983l) . iHevvaerts fc Priest! l|l984ri 
and others. ICargilll l|l994|) pointed out that impulsively 

heated loops would cool by conduction and radiation with 
observable results. The general view of these studies was 
confined to the problem of releasing magnetic energy in 
the corona with the goal to keep it hot. 

This is not supported by recent observations. Deep ex- 
posures of the soft X-ray emission in a quiet region by 
Yohkoh/SXT revealed that the emission measure of the 
corona in certain pixels was not constant over t he ob - 
serving time of about one hour l|Krucker et all Il997f) . 
It suggests that the material content of the low corona, 
where most of the soft X-ray emission originates, varies. 
More sensitive observations of EUV lines from the corona 
show variability of the majority of pixels (85%), in cluding 
some in the faint intra-ce ll regions ( Benz fc Krucker. 1998: 
iBerghmans et all Il998(l . The increases of the emission 
measure are not adiabatic compressions and can only be 
interpreted as additions of new material into the corona. 
The heated material subsequently seems to cool on a time 
scale of the order of 15 minutes. Thus the coronal ma- 
terial in the lower corona appears to be not heated, but 
rather continuously replaced. We may add here, however, 
that a heating process could still be hidden in the quasi- 
stationary background of the emission measure. 
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Some large ev ents have recently been tested by 
iBrown et all i200(^ for their relation between increases 
in emission measure and temperature. The authors con- 
clude that this relation corresponds more to an impul- 
sive heating of the upper chromosphere and subsequent 
expansion rather than a heating at the top of a coro- 
nal loop followed by conductive readjustment. Thus the 
observed increases of the coronal emission measure may 
be interpreted as " chromospheri c evaporation" similar to 
reg ular flares in active reg ions. iBenz fc Kruckeil l|l999t) 
and lKrucker &: Beii3 l)200Cl|) find very little difference be- 
tween flares and the observable properties of the emis- 
sion measure increases, and thus confirm the latter as real 
nanoflares. In the f ollowing, we w ill use the term nanoflare, 
first introduced bv IParkenl) 1 983(l as a theoretical concept, 
for any brightening of the quiet corona with energy below 
ap proximately 10^^ erg. 

iKrucker fc Bend l)l998l) estimated the total energy in- 
put by the emission measure increases to be 16% of the 
calculated total radiative loss of the observed region. This 
number depends strongly on the sensitivity of the instru- 
ment and some model parameters. In particular, the ef- 
fective line-of-sight thickness of the coronal plasma (or 
height for observations in the center of the disk) can- 
not be measured and must be assumed. The distribu- 
tion of the events in energy is therefore still contro- 
versial. Most observers report a power law shape, but 
widely disagree on the expon ent, which ranges from —1.45 
ijBerghmans fc CletteL f 19991 measuring radiation loss in 
Fe Xin to -2.59 dKrucker fc B enz. 1998, measuring emis- 
sion measure increases V iBenz fc Krucker (,2001a.') reported 
agreement between their EI T based analysis and stud- 
ies based on TRACE dat a bv lParnell fc JupdI |)2000|) and 
lAschwanden et ahl l)200(]|) . if the same method is used. For 
a flare model with a height proportional to the square root 
of the flare simultaneous peak time within 2 min- 

utes over the flare area and no further flare selection, all 
three investigations yield a power law index in the range 
between —2.0 and —2.4, the most likely new EIT value 
being —2.3. Clearly, nanoflare heating of the quiet corona 
is strongly supported by the above observations on total 
energy input and the relatively steep slope of the energy 
distribution, which suggests that the smallest flares con- 
tribute most to the heating. 

Here we address the question of whether nanoflares 
can account for all of the observed properties of the emis- 
sion measure in the quiet corona, including (i) the gen- 
eral appearance of individual pixels' emission measure 
time dependence, including background and nanoflares, 
(ii) the absolute value of the quasi-steady emission mea- 
sure (equivalent to the total radiative loss), and (Hi) the 
average Fourier spectrum of the emission measure in time, 
reported to be a power law with an exponent of —1.76 
l)Benz fc Kruckeilll998|) . 

First, we derive the conversion of flare frequency dis- 
tributions in energy for different assumptions, in partic- 
ular the conversion between distributions of flares, where 
pixels have been grouped into events, and distributions 



of single pixels, as well as the conversion between differ- 
ent height models (Sect. In Sect. |2| we simulate the 
emission measure fluctuations in the quiet corona by tak- 
ing the observed flare frequency distribution in energy, 
making certain assumptions on the spatial and temporal 
shape of the events and extrapolating them to the small- 
est energies needed to explain the total radiation loss from 
the corona. By choosing an adequate time proflle for each 
flare, we can simulate the time dependence of the emission 
measure for an arbitrary pixel. We also obtain an averaged 
power spectrum. In Sect. ^ we calculate the expected value 
of the power spectrum analytically and in Sect. |31 the re- 
sults are discussed. Section El summarizes the conditions 
needed for the nanoflare model to reproduce the observed 
features of the variations and quasi-steady background of 
the quiet corona. 

2. Flare Distributions 

The observations that have been made of the quiet Sun 
to determine the energy input into the corona are based 
on a grid of many pixels in a large field of view. For 
each pixel and time interval , the emission measure is esti - 
mated (cf . iBenz fc Kruckert Il998t IParnell fc JuppL l2000|) . 
Neighboring and simultaneously occurring events are com- 
bined and interpreted as one flare. The thermal en- 
ergy of an em ission measure increase can be estimated 
l|Benz fc Krucker. 1998,) 

E « SkeTVAMAh, (1) 

where the emission measure increase AM, the flare area 
A and the temperature T are observed quantities and fcs 
is the Boltzmann constant. Equation assumes that the 
newly injected material is much dens er than the pre-flar e 
density in the loop. As pointed out bv lBrown et alJl|200nl) . 
this assumption yields an upper limit to the energy. On 
the other hand, the flare energy release may also directly 
heat the coronal material, whose energy is not included in 
Eq. QJ. In the following, we will use Eq. as a rough 
estimate for the thermal energy input by a nanoflare. The 
height (thickness) h of the flare loop lies in the line of 
sight and is therefore not easy to determine. We will use 
it here as a free model parameter that also includes the 
filling factor. Thus h denotes the effective height. 

All observers agree that the flare frequency distribu- 
tion function depends on the energy as a power law 

f{E) = foE-', (2) 

where the flare distribution is defined as the number of 
flare events within a given energy interval divided by the 
total observed area, total observation time and energy in- 
terval. Similarly, t he flare distribution i n area was reported 
to be a power law ((Aschwanden et al.[EoO(ilAlietti et all . 
1200(1) . Supported by these observed distributions, it is as- 
sumed in the following that flares are self-similar in ther- 
mal energy. Thus on average, a fiare is completely charac- 
terized by its energy. This assumption will later be tested 
against observations. 
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Let us therefore write the emission measure increase, 
the flare area and the flare height, each averaged for a 
given thermal energy, in terms of some power of the energy 



AM = c-E'' 
A = a-E" 
h = b-E^. 



(3) 
(4) 
(5) 



Where no confusion is possible, we will use the same sym- 
bols for averaged entities, such as AA^, and measured 
quantities in individual events. Equations ®-© together 
with Eq. 1^ imply the two conditions 



a + /3 + 7 = 2 
abciSkeTf = 1. 



(6) 
(7) 



It is seen immediately from Eq. that the energy 
of a flare depends on the height assumption, and so do 
the parameters /o and i5 in Eq. ij^J. In the first and sec- 
ond part of this section we study the scaling relations and 
conversions between flare energy distributions for different 
assumptions for h as well as for time variations of individ- 
ual pixels. In the third part of this Section we will apply 
these relations to numerical values from observations. 

2.1. Distribution in flare area and per pixel 

Because both the emission measure increase and the area 
of a flare are observable quantities, we write 



AMiA) = mA''. 



In terms of the above variables we have 

7 
a 



m 



c ■ a 



(8) 

(9) 
(10) 



There are four equations (Eqs. IHl El El ^^nd I10|l and six 

variables (a, b, c, a, /3, 7), thus the system is still not 
fully determined. This is because we have the freedom to 
make an assumption regarding the height dependence. A 
model for the flare height fixes b and (3, and the system 
of equations can be solved. There is one more difficulty, 
however, and that is to determine /i and m. As Eq. (jH]) is a 
direct relation between variables of individual events and 
not a distribution, the scatter is large, which is obvious 
in Fig. ^ There is another way to determine /i and to, 
using the comparison between flare distribution and the 
energy distribution of emission measure increases observed 
in individual pixels. In the following, we refer to the latter 
as pixel distribution and mark it with the index p. 

The pixel distribution function accounts for events in 
pixels without adjacent pixels being grouped together to 
form one flare. Thus for flares covering more than one 
pixel, the pixel distribution is the flare distribution times 
the average flare area measured in units of pixel area Ap 



U{Ep)dEp = ^f(E)dE, 



(11) 



where the thermal energy of a pixel is Ep — 
'iksT AM.pAph. The average emission measure increase 
in a pixel is simply AAip = Ap/A- AM.. It follows that 

Ar. 



E„ 



A{E) 



E. 



(12) 



We insert Eq. Q into Eq. (fT^ and solve for E to obtain 
an average relation between the thermal energy input in 
flares and in pixels 



E = 



a 

Jip 



(13) 



Taking the derivative 



1 



\~a\A 



^ eT- 



dE 

rewriting Eq. Hll|) as 

fpiEp)dEp = ^^^mE,])§-^ 

and inserting Eqs. (0), (jH), lO and lO into Eq. (flS)) . 
leads to the flnal result for the relation between the flare 
distribution and the pixel distribution 

2-5 



(14) 



(15) 



fp{Ep) - 



1 



1 — a 



.foEp 



(16) 



The exponent a and the factor a then are 



a 



Arn 



fpO 



Sp-2 



(17) 



(18) 



Equation ^ together with the height assumption used 
to derive the thermal energy then yields 7, and /x follows 
immediately 

- m - 2) 



7 Sp + S-^ 

= -= 7 

a Oj 

Similarly, one derives 



S 



biSkeTy 



(19) 



(20) 



2.2. Scaling relationships 

The primary free model parameter is the height of the 
flare (Eq.[SJ), i. e. the parameters b and /3. They determine 
the thermal energy by Eq. (© when T, AM and A are 
observed. In the following we derive the conversion be- 
tween the flare frequency distributions of different height 
dependencies. 

Let the new flare distribution be /'(£"), where E' = 
ZksT^ AMAh' . Then the conversion from the old flare 
distribution (Eq. \^ to the new one is given by 



nE')dE' ^f{E[E'])^dE'. 



(21) 
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It is seen immediately that the identity El\fh — E' /^/h! 
holds. Inserting Eq. jSJ, we obtain from it the condition 



Solving for E leads to 

/ , 2-13' 

E{E')^^-j E'—. 

The derivative is 

dE ^2~(3' fb\^ ^^ 
dE' 2-/3 \b' J 



(22) 



(23) 



(24) 



Inserting now Eqs. if^ and (Plj) into Eq. lf?T|l gives 
the new flare distribution 



„ 2-3' fb\^ I 6(2-l3')-f3 + l3' 

f'iE') = ^ (- j foE ^ 
defining the new power law exponent S' and /q. 



(25) 



2.3. Numerical values 

The observed quantities are t he emission measure, the 
flare area and the temperature. iBenz fc Kruckej l)2001b(l 
flnd an average temperature of 1.46 • 10^ K for the 23 
largest nanoflares that occurred in a quiet region within 
one hour. The temperature does not depend on the flare 
energy, and its distribution is narrow. The highest ob- 
served value was reported to be T = 1.63 • 10^ K. We use 
the former value for our models and will discuss it more in 
Sect. El The pixel area is given by Ap — 1900 kmx 1900 km. 

Making an assumption on the flare height, the flare 
energy, flare distribution, pixel energy and pixel distribu- 
tion can be reduced, and thus also the parameters m and 
H defined in Eq. JSjl according to Sect. 12. ll Once we know 
m and /i, we can solve for the parameters a, b, c and a, 
/?, 7. Also, we can calculate a new flare distribution from 
Eq. H25|) and a new pixel distribution from Eq. I|16l) for a 
different flare height. This shows the effect of the model 
assumption on the distribution and characteristic ener- 
gies. Finally, the biggest interest is in the parameters for 
a model with flare height that scales with flare size. 

Table n displays some parameters of interest. The val- 
ues for /o and 5 for h = 5000 km are from lKrucker &: Bend 
l)l998|) and fpQ and Sp for h — 500 km are from Krucker 
(private communication) . All the other values from /g to 7 
are calculated using the methods of the previous sections. 

First, /o (500 km) and 5(500 km) have been calculated 
from Eq. if^ . then m and /i from Eqs. and ifT^ . 

The relation of Eq. ||SJ| is displayed in Fig. ^ (solid line). 
T his figure also contains the largest 22 events reported 
bv lKrucker fc Benz ( 2000 ) from a large field of view, thus 
the emission measures are relatively large (their Table I, 
shown here by -I-), and the largest 6 events from a smaller 
set (their Table II, shown by *). The calculated relation 
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Fig. 1. Emission measure per flare versus flare area. The 
crosses are the largest flares observed in a large field of 
view and the asterisks refer to the larg est events in a small 
field of view iKrucker & BenA\200(l} . The solid line rep- 
resents the relation AM = mA^^ (Eq. using the pa- 
rameters calculated from Eqs. J^19\) and Y20\) for (3 = 0. 



between emission measure and flare area matches well that 
observed for individual events. Knowing the parameters 
m = 7 • 10^° and /i = 1.34, which are independent of the 
height model, all the other parameters can then imme- 
diately be calculated for any chosen height dependence. 

Rather large nanoflares are seen to cover an area of 
approximately rectangular shape with a length to width 
ratio q = l/w of about 10. Thus 



A = l' 



(26) 



As ordinary flares are believed to be generated in loop- 
shaped magnetic fleld structures, nanoflares may be ex- 
pected to originate in similar geometries. Thus, we as- 
sume the nanoflares to occur in torus-shaped loops with 
a loop-thickness equal to the flare width. In our approx- 
imation this loop-thickness shall just be the unknown ef- 
fective height, thus h = w = ^ (A/10) and (3 = a/2 and 

& = 7(^710). 

The largest single-peaked flare observed was derived 
to have an energy of i^maa; (5000 km) = 2 • 10^^ erg 



llKrucker fc Bend. 120001 Table I). The minimum energy 
can be derived by equating the total radiation output 
Plot — 4.5 • 10^ erg~^ cm~^ s~^ with th e integrated flare 
energy input. This value was derived bv lKrucker fc Bend 
(1998) from the total emission measure observed in coro- 
nal EUV lines, averaged over the whole field of view, and 
includes all radiation losses in the continuum and the lines 
from UV to X-rays. Solving 



Pf.nt — 



E ■ f{E)dE 



(27) 
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Table 1. Observed and derived parameters of nanoBare 
energy distributions for different height models given by 
h, b and (3. The temperature used was T = 1.56 • lO^K. 



h 


5000 km 




500 km 


^/(A/10) 


fa 








3.38 ■ W 


5 


2.59 




2.59 


2.31 


J 


■|^gl06.12 




2^q103.6 


1.30 ■ 10^° 


On 


6.04 




6.04 


3.05 


a 


1.8 ■ 10"" 




4.7 ■ lO^'' 


0.209 


b 


5 • 10* 




5 • 10^ 


0.145 


c 


2.7- 10" 




1.0 • lO^'^ 


7.91 • 10^^ 


a 


0.854 




0.854 


0.704 


p 










0.352 


7 


1.146 




1.146 


0.944 


Emin 


2.2 • 10"'' erg 


9.8 


• 10"' erg 


3.7 • 10""' erg 


Eres 


9.0 • 10^"* erg 


2.8 


■ 10^* erg 


3.1 ■ 10^'' erg 


Emax 


2.0 ■ lO^'^ erg 


6.3 


■ 10^^ erg 


1.3 ■ 10^^ erg 


Imin 


1230 km 




530 km 


250 km 



one obtains for 5 > 2 

Emin ~ ( — 7 Ptot] " ■ (28) 

V /o / 

Here we have used the fact that for our parameter range 
the upper boundary term of Eq. I|27|) is vanishingly smah. 
The length of the smallest flares is given by Eqs. I@J and 

- (10 • a • • (29) 

If the flare size is equal to the area of a pixel, Eq. (|12|l 
shows that the flare energy is equal to the energy in 
the pixel increase. At this energy, Eres, the flare energy 
distribution will have a low-energy roll-over due to sub- 
resolution events. The energy Eres is obtained by solving 
Eq. H13|) for E = Ep = Eres- It is co nsistent with the 
cutoff energy of 9 • 10^'' erg observed by iKrucker fc Bend 
l)l998l Fig. 2, assuming h = 5000 km). 

In Fig. [5] the different flare distributions as well as the 
pixel distribution for h — ^ (A/IO) are shown for com- 
parison. The distributions are displayed for energies be- 
tween Eres and Emax sls they were observed. The pixel 
distribution only makes sense for energies higher than the 
resolution energy, because only these flares cover one or 
more pixels. This can be illustrated by the following con- 
sideration: Let us assume a flare with energy smaller than 
Eres, thus the flare area is smaller than the pixel area. 
According to Eq. Hll(l . the number of pixel events in the 
energy interval dEp would be smaller than the number 
of flare-events in the corresponding energy interval dE. 
This contadicts the assumption that, when any such flare 
occurs, it has to appear in exactly one pixel. We neglect 
here the effect of subresolution flares spreading over pixel 
boundaries. Thus, for flares with area/energy smaller than 
the resolution, the pixel distribution is equal to the flare 
distribution and not given by Eq. (|11|) . which in that case 
makes no sense. This can also be visualized by the defini- 
tion for the flare distribution, which is the number of flares 



pixel 




Fig. 2. Various flare distributions and the pixel distribu- 
tion for h ^ ^/A/T0 are shown for different height as- 
sumptions. The bold line labeled with h=5000km cor- 
responds to the flare energy distribution observed by 
\Krucker & Bend {1998). The other distributions are de- 
rived using Eqs. ^J6]) and i)25|) . The asterisks denote Eres 
and Emax. 

divided by total observed area, time and energy interval. 
Obviously, for flares within the observed area this fraction 
is constant. Reducing the area to the size of a pixel doesn't 
change the flare frequency distribution for flares with still 
smaller size. 

3. Simulation 

In this section we simulate numerically the emission mea- 
sure variations in individual pixels. The goal of the simu- 
lation is to find the conditions for making the simulated 
pixels resemble the observed emission measure in form 
and background level, and to make the averaged simu- 
lated power spectrum match the observed one in slope 
and absolute value. 

It is assumed that the total emission measure in a pixel 
at a certain time is given by the sum of the emission mea- 
sures of all flares brightening at that time. The flares are 
described by their energy, which is distributed according 
to the observed pixel distribution for flare energies in the 
observed range Eres to Emax and flare distribution for 
the extrapolated sub- resolution events from Emm to Eres ■ 
The emission measure increase is given by Eq. Because 
this increase was defined as the difference between a max- 
imum and the preceeding minimum of a brightening, it is 
about equal to the maximum value of the time dependent 
emission measure of a flare. The emission measure of a 
flare increases up to this value and then decreases again 
according to the flare's time profile. 

In the simulation, the flares are randomly distributed 
in time, according to the flare and pixel distribution. Each 
flare is defined by its maximum value and the time pro- 
file. For every observation, the total emission measure then 
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is the sum of the respective flares' emission measures at 
that time. To make adequate comparison with the obser- 
vations, we choose the total observing time tabs and the 
equally distributed number of observin g points Nnh.-< in 
the same way as in the observations bv iBenz fc Kruckeil 



l|1998(l . Of the total time dependent emission measure, the 
power spectrum is taken numerically. This process, simu- 
lating a pixel and taking the power spectrum, is repeated 
many times and the resulting power spectra are averaged. 
The results, the single pixel emission measure level and 
variations and the average power spectrum, are then com- 
pared with the observations. 

3.1. Time profile 

The power spectrum of the emission measure variations 
in time per pixel depends on the time profile of individual 
events. Let Ai be the emission measure of one flare. The 
assumed self-similarity of flares requires that M must de- 
pend on the average only on energy in amplitude as well 
as in duration. The temporal evolution of a flare is given 
by the rate at which plasma is heated to coronal temper- 
atures, expands and cools. We will assume that the dura- 
tion is proportional to E'^ . The average emission measure 
of flares with energy E is therefore 



M{E,t) = cE''g{E-H), 



(30) 



where g is the time profile. Note that the maximum value 
of the emission measure of one large flare is assumed to 
be the observed emission measure increase. Therefore the 
time profile g{x) = g{E~'^t) must be chosen in such a way, 
that max{g) — 1. The emission measure increase IS.M — 
max{Ai) — cE^ then follows the relationship Eq. 

Let t{E) be the characteristic duration of a flare with 
thermal energy E. We then have 



t{E) (x^ ocE". 
E 

The emission measure M is proportional to the observed 
energy loss rate from radiation. The corona may also lose 
energy from conduction at the foot-points of the loop with 
a rate Q. Then E oc Ai + Q. Assuming that Q follows the 
same power law as M or is much smaller and therefore 
can be neglected, we obtain 



E 



(32) 



AMiE) 

Inserting Eq. (|2Jl and comparing with Eq. H31|) , we imme- 
diately obtain 5 = 1 — 7. The simulation results bear out 
this assumption. 

3.2. Density 

In the derivation of the equation for the thermal energy 
of a flare (Eq. ^ it was assumed that the density n{E) is 
approximately constant in the flare volume. Thus, 



where V — A ■ h = abE°'^^. Solving for n and using 
Eq. lO, we get 

"-(^)^"^- 

For 7 < 1 (s > 0), the density decreases with energy 
as the duration increases. A physical explanation may be 
that the longer life time of large flares is partially a result 
of lower density. Note that the effect is small as s is close 
to zero in all mode l s. A s imilar result has been derived by 
lAschwanden et al.lll200nl) . 



3.3. Results 

In the simulation we use the area dependent height model 
h = \J A/10 with the parameters derived from observa- 
tions (Table ^ S"^"^ column). We just remark here that 
this is not the only height assumption which yields rea- 
sonably good results. In particular, the observations are 
also reproducible for a model with constant height with 
only minor adjustments in other parameters. Note that 
the factor c, which is directly proportional to the emission 
measure (Eq. and therefore proportional to the square 
root of the power spectrum, is inversely proportional to 
the temperature squared. Thus, the absolute value of the 
power spectrum is proportional to the minus fourth power 
of the temperature. For fine tuning the absolute level of 
the power spectrum we will adjust the temperature in 
the simulation. This is equivalent to saying that the ef- 
fective thermal energy is not given by Eq. but rather 
by = eSkBTobsV AMAh, with Tabs the observed tem- 
perature and e a correction variable for the true thermal 
energy. Thus, the temperature we use is T = eTobs with 
Tobs = 1.46-106K. 



The following time profile is used 



(31) 9{x) 



exp(^a:), a; < 
exp(— 77a;), a; > 



(35) 



AM := / 
Jv 



dV 



(33) 



with X = ti the peaking time of the i-th flare 

with energy E and duration r = t{E) = {E/EoYto, 
tq — 1000 s the characteristic time of a flare with energy 
Eq = 5.8 • 10^^ erg, 5 = 1 — 7 = 0.0556 and ^ and 77 pos- 
itive free parameters. For the presented results we have 
used ^ = 20 and rj = 0.8. The emission measure per pixel 
is defined as EM := M/Ap and is of order of 10^^ cm~^. 
The pixel area is Ap = 3.61 • 10^^ cm^, and the total obser- 
vation time is tobs — Nobs At with Nobs = 21 the number 
of observations and At — 121.9 s the time step. 

In Fig. 131 some simulated pixels are displayed. The dot- 
ted curve is the summed emission measure over the flares 
in the sub-resolution regime. We see that they are respon- 
sible for the background level. The dashed curve occurring 
occasionally is a resolved nanoflare. The solid curve is the 
sum of the two: the total emission measure per pixel. 

To simulate the effect of noise, the total emission mea- 
sure at an observed time is distributed randomly within 
the emission measure noise interval clEM . Figure 0] shows 
the effect of noise for the same model parameters. 
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Fig. 3. Emission measure per pixel (solid curve) displayed 
for random pixels. The dotted curve is background from 
sub-resolution flares. The spatially resolved nanohares are 
shown by a dashed curve. 



Figure shows the averaged power spectrum of pixels, 
simulated with the same parame ters as in Figs . 01 and H] 
The solid curve was observed by iBenz fc Krucker 1,1998.) 
averaged over all pixels, and here corrected in the ordinate. 
The slope of the simulated power spectrum (dotted curve) 
at low frequencies is mainly given by the choice of the free 
parameters ^ and rj (see also Sect.EJ. The noise increases 
the power at all frequencies (dashed-dotted curve) . At low 
frequencies the slope is not changed, but at high frequen- 
cies the curve is considerably bent upward, consistent with 
observations. The absolute values of the power spectrum 
depend strongly on the temperature as discussed in the 
beginning of this subsection. The best fit is found to be 
T « 1.56 • 10*^ K, thus e « 1.068. 

The averaged power spectrum was found to be highly 
sensitive to the time profile of a flare, which is poorly 
known from direct observation. 



4. Analytic power spectrum 

In this section the power spectrum of a given nanoflare 
model distribution is calculated analytically. 

Let the expected value of the averagedj)Ower spectrum 
at frequency oj be £{\Mtot{(^)\'^), where Mtot denotes the 
Fourier transform of the emission measure of a pixel. If the 
flares follow a Poissonian distribution, are independent of 
each other and identical at the same energy, it can be 
shown that the total power spectrum is composed of a 
flare component and a noise component (see Appendix IXjl 

£{\Mtot{uj)\^) = £{\Mfl{i^)\^) + £{\Mno^se{i^)\^) (36) 

with 



£{\MfiiLu)\')=Ap fiE)\M{E,Lu)\'dE (37) 



Fig. 4. Noise added emission measure per pixel (solid 
curve) displayed for random pixels (different from Fig.\^. 
The dashed curve is the emission measure without noise. 



and 



At 



(38) 



dA4 being the noise level. Note that in the simulation as 
well as in the observation the numerical Fourier transform 
is normalized by division of the total observing time tabs-: 
so here we use the same normalization. Again we consider 
pixels. Thus the pixel distribution has to be taken for large 
events, and the integral splits into two parts. 



Art 



f{E)\M{E,Lu)\^dE+ 



fp(Ep)\Mp{E,Lo)\^dEp 



(39) 



With fpdEp = A/Ap ■ fdE and Mp = Ap/A ■ M we get 



£(\Mfi(co)\') = Ap 



f{E)\M{E,Lu)\'dE+ 



^f(E)\M{E,u:)\'dE 



lE,^.. A{E) 

The Fourier transform of a single flare is given by 

M{E,uj) = cE^gitlr) 
= cE'^Tg{TUj). 

With g{x) from Eq. (|35() we have 
.9(fc) - 



(40) 



(41) 
(42) 



(43) 



(e-ifc)(77 + iA:) 
with k = TLO and 

\M(E,uj)\'' ^c^E''^T{Ef— , , , , 

(C^ + (TCj)2)(7y2 + (T-t^)2) 

Using uj = ^-nv and inserting all into Eq. 1)40(1. we obtain 



(44) 
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frequency [mHz] 



1.0 

frequency [mHz] 



Fig. 5. Simulated power spectrum averaged over 10000 
pixels with ^ = 20, r] = 0.8. The dotted curve is the power 
spectrum without noise and the dashed-dotted curve the 
one with noise, the latter is to be compared with the ob- 
served power spectrum (solid curve). 

Si\MM,.r)^A,c'^M^ + vr- (45) 

^£ (e + (2^T„ (f )%)2)(r;2 + (27rro (f ^ 

Ap ^"F E'^-^-"dE 

which can be computed numerically. The solution is shown 
in Fig. 1^ dashed curve, with all parameters the same as 
in Sect. 13.31 in particular ^ = 20 and 77 = 0.8. In the dot- 
ted curve, aliasing is additionally taken into account. The 
noise term £ {\A4noise{^)\'^) is added to the aliased power 
spectrum (dashed-dotted curve) with dAi = Ap ■ dEM. 
The solid curve shows the observed power spectrum. 

5. Discussion 

It is clear from the previous two sections that a model 
exists that can explain the emission measure and its fluc- 
tuations observed in the quiet solar corona. In this section 
we address the question of how sensitive the results are 
to the choice of parameters. Assumptions were made on 
the height dependence of the thermal energy, extrapola- 
tion of the flare distribution to energies several orders of 
magnitude smaller than the observed range and the shape 
and energy dependence of the time profile of single pix- 
els. Also, an assumption used throughout this model is 
the self-similarity of flares at different energies in the av- 
erage, which is further supported by having found a set 
of parameters that explains all the available observations. 
In the following considerations, we will focus on a model 



Fig. 6. Analytically calculated power spectrum with ^ = 
20, 77 = 0.8. The dashed curve shows the true power spec- 
trum, the dotted curve the power spectrum with alias- 
ing and the dashed-dotted curve the power spectrum with 
aliasing and noise. The solid curve is the observed power 
spectrum. 



with h — A/10, which seems to be a plausible choice. 
However, it should be mentioned that this model is not 
unique in meeting the above requirements. In particular, 
models with a constant height of 500 km or 5000 km, re- 
spectively, also yield reasonably good results. The remain- 
ing free parameters for fitting ,T,^, and rj, however remain 
practically the same. 

The three main observational constraints mentioned in 
the introduction are the typical emission measure fluctu- 
ation of a pixel, its absolute value in emission measure 
and the shape and value of the power spectrum. As seen 
in Fig. 121 the sub-resolution flares contribute most to the 
background level of the emission measure. If we raise the 
minimum energy or, equivalently, reduce the total input 
power, the required number of small flares decreases and 
the background emission level drops drastically. To ex- 
plain the observed level in emission measure we need to 
extrapolate the energies to Emin- The background level 
of the preferred model yields values at the low end of the 
observed emission measure in the quiet corona. A con- 
stant emission measure background (e. g. produced by an- 
other heating process) could be added without changing 
the power spectrum. 

The nanoflare time profile plays an important role in 
all three aspects. The background level in emission mea- 
sure depends on the duration of a flare. A slow increase 
and decrease in the time profile makes a flare contribute 
longer to the background, whose level then rises. In this 
case the fiuctuations become smaller. For a fast increase 
and decrease the background level drops and the fluctua- 
tions become more peaked and accentuated. 

The rate of change has a strong effect on the slope in 
the power spectrum. For the time profile given by Eq. (|35|l . 
the slope of the power spectrum is mainly determined by 
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frequency [mHz] 



Fig. 7. Analytic power spectrum for a symmetric 
nanoRare time proEle with ^ = rj = 3.41. The dashed 
curve displays the basic power spectrum, the dotted curve 
the power spectrum with aliasing and the dashed-dotted 
curve the one with noise and with aliasing. This final cal- 
culated curve drops drastically at higher frequencies, not 
even noise can lift it enough to match the observed spec- 
trum (solid curve) at high frequencies. 



the free parameters ^ and 77. A slope with the observed 
value K = —1.76 between the lowest two frequency points 
can be obtained by a wide parameter range, e. g. 77) = 
{(3.41, 3.41), (10, 1.2), (15, 1.0), (20, 0.8), (50, 0.7)}. For a 
symmetric time profile, i. e. ^ = 77 = 3.41, the power 
at high frequencies drops so drastically, that a bump ap- 
pears at medium frequencies that cannot be flattened out 
by noise (Fig. [7Il. The more asymmetric the time profile, 
the straighter is the slope of the power spectrum and the 
more it resembles the observed one. That the time profile 
should be so asymmetric is surprising, although it is con- 
firmed immediately by Eq. (|45|l for the time profile used. 
If ^ is large compared to the dependence, one bracket- 
term in the denominator of the integral is almost negli- 
gible. If, additionally, rj is small, the power spectrum is 
approximately proportional to iy~^. Added noise reduces 
the steepness of the slope further which then matches the 
observed k. The small discrepancies in the level of the 
simulated and expected power spectrum come from its 
sensitive dependence on the parameters a, b and c and 
numerical errors in the last digits thereof. A small depen- 
dence of K also comes from the maximum energy Emax- 
The larger the maximum energy, the steeper the slope. 

6. Conclusion 

Deriving distributions of fiare energy and pixel variations 
from observations requires some assumptions, in particu- 
lar a model of the effective line-of-sight thickness (height). 
We have here developed the analytical tools for the trans- 
formation between different assumptions published in the 
literature. These tools are the basis on which we develop 



a numerical model that simulates the time behavior of 
individual pixels and nanofla res. It is used here to r epro- 
duce the FIT ob s ervati ons of lBenz fc Kriickerl ((l998l) and 
iKrucker fc ReuTl l)l 99Sl) . 

We have approached the question of nanoflare heating 
by searching for a model of energy input by nanoflares 
that can reproduce the three constraints given by the ob- 
servations of individual pixels: the general appearance of 
the emission measure variation in time per pixel, the av- 
erage level of the radiation loss, and the average Fourier 
power spectrum in time. 

We obtain the power spectrum in two different ways, 
averaging the simulation and calculating the expected 
value directly. As the evaluation of the integral con- 
sumes much less calculation time than the simulation, this 
method allows us to test systematically different time pro- 
files. 

The analysis by analytical and numerical studies has 
shown that there is at least one such model, assuming 
self-similarity in energy. Closest to reality may be a model 
with a height proportional to the square root of the area. 
The model parameters correspond to the last column in 
Tabled To reproduce the observed radiation loss, the ob- 
served range of nanoflares (having a lower limit at about 
3 • lO^'' erg, cf. Fig.lJJ needs to be extrapolated to lower en- 
ergies by four orders of magnitude (Table T|). Th is energy is 
far below the model suggested bv lParke" l)l983() . although 
there is no stringent theoretical lower limit of magnetic 
energy release. Nevertheless, the extrapolation is hypo- 
thetical for several reasons, in particular the assumption 
of self-similarity implying a constant power-law index for 
the distributions in energy, area and height. Other height 
assumptions, e. g. constant height, also yield reasonable 
results after slight adjustments for some parameters. 

The minimum energy depends on several assumptions. 
It is larger if the fiares introduce other forms of energy 
into the corona than just the thermal energy of evapo- 
rated material, such as fiuid motions or wave energy. They 
would heat the material already existing in the corona 
and yield a constant background emission measure. Such 
a background level would reduce the thermal energy in- 
put requirements by nanoflares and enhance the minimum 
energy. A test on the importance of a steady emission 
measure versus flares could be made by observations of 
much higher sensitivity than possible today with EIT and 
TRACF. 

A hint of a break-down of our self-similar model may 
be the small length of fiare loops derived (Table A 
possible remedy may be to assume an energy dependence 
of the length to width ratio q. One may expect from a 
reconnection scenario that the width and effective height 
decrease faster with energy than the length, thus forming 
the smallest flares in thin, long loops. There are currently 
no observations of this ratio over a sufficient range of flare 
energies. High resolution observations could clarify this 
point in the future. 

The fitting model is not unique, particularly in the 
choice of height assumption and average time profile of 
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nanoflares. Nevertheless, we found that the three obser- 
vational constraints severely limit the range of free pa- 
rameters once the model assumptions have been made. In 
particular, the power spectrum is found to be very sensi- 
tive to the chosen time profile. As the exact shape of the 
time profile is not easily observable, this can be used to 
test how nanoflares evolve and disappear. 
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Appendix A: Derivation of the expected value of 
the power spectrum 

Consider the observation time interval tohs divided into 
subintervals of equal length through the points to — 
0,ti, . . . ,t]sf — tabs, with N ^ Nobs- If the subinter- 
vals are chosen to be sufficiently small, we can assume 
without loss of generality that each flare peaks at one 
of the points t„. Further let us divide our energy range 
into M equal intervals of length AE through the points 
Eo = Emin, El, . . . , Em = Emax- Let pmn dcnotc the 
number of flares of energy between Em and Em + AE 
peaking at t„. For any to, let Pmri be Poissonian, indepen- 
dent and identically distributed. Their variance is then 
equal to their expected value and we have 



Var{pmn) = £{Pnin) 



f{Em)AE 

Jf{E)dE 



(A.l) 



with f{Em) the frequency distribution of the flares, nor- 
malized to 



f{E)dE 



N 



AptQ. 



p^obs 



(A.2) 



Let the emission measure of a flare with energy Em 
peaking at time i„ be 

Mmnit)^M{Em,t-tn). (A.3) 

The total emission measure at time t is then given by 



Mfiit) 



m.n 



Pm7i-^^m7i (^) ■ 



Thus, we have 



Mfiiiv) 



m.n 



therefore 

\Mfi{uj)\^ = M-M* 



(A.4) 



(A.5) 



(A.6) 



— ^ ^ Pm 1 n 1 Pm2 n2 

Tl]^ ,712 

M{Em, , uj)M* {Em, , t^)e-''^(*"i ) . (A.7) 

As the expected value operator commutes with the sum, 
we have 

£i\MMio)\') = 

^ ^ ^{Pm\n\Prti2n2^ 

mi ,m2 
ni ,712 

M{Em, , u;)M* {Em, , uj)e~''^^'-^ ) (A.8) 

mini ) + £{p mini m2n2 )\ 



ai ,7712 

■11 ,T12 



M{Em, , uj)M* {Em, , c^)e-''^(*"i ) . (A.9) 



Inserting Eq. (|A.ip 
S{\MfM?) = 



■Aptohs 

N 



J2 f{Em)M{Em, Lu)M*{Em,u;)AE + 



+ f{Em,)f{Em2)M{Em,,Uj)M*{Em2,Uj)AE^ 



{Aptob 



iV2 



E 



-iw(t„j-t„2) 



(A.IO) 



and taking the sum over ni, n2 and n, respectively, one 
obtains 

£{\MMi.r) = 



Aptobs J2 f{Em)\M{Em, uj)\^AE - 
m 



N 



(A.ll) 



f{Em, )f{Em2 )M{Em, , io)M* {Em, , Uj)AE^ 



For N oo the second term vanishes. Taking the limit 
AE ^ we finally obtain 



£{\Mfi{u;)\') Aptobs / f{E)\M{E,cj)\'dE. (A.12) 

White noise adds in the following way. Let the total 
emission measure with noise be 



Mtot{t) =Mfl{t)+Mnotse{t) 

where 

Mnoise{t) = rdM 



(A.13) 
(A.14) 
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with r a random variable between —0.5 and +0.5 and dM Ulmschneider, P., Rosner, R., & Priest, E. R. 1991, 

Mechanisms of chromospheric and coronal heating 
(Springer- Verlag, Berlin) 



the noise interval. Then 

Mtot{i^) = Mfi{Lj) + rdM 
\Mtotf = \Mfif + r^\dMf + 

+r (MfidM + M}idM) 

£{\Mtot?) = £(\Mfi\')+£{r^)\dM\^ + 

+£{r) (MfidM* + M}idM^ . 

It is £{r) = 0, £{r^) = 1/12 and 



Nobs 

= tobsAt{dMf. 

Thus, 

^(l^totHl') = 



(A.15) 
(A.16) 

(A.17) 

(A.18) 

(A.19) 
(A.20) 



^ohs 



/At 
f{E)\M{E,oj)\^dE + —{dM)' 



(A.21) 
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